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ABSTRACT 

We present adaptive optics photometry and spectra in the JiJ i^i-bands along with high spectral 
resolution X-band spectroscopy for each component of the Z Canis Majoris system. Our high angular 
resolution photometry of this very young (<1 Myr) binary, comprised of an FU Ori object and a Herbig 
Ae/Be star, were gathered shortly after the 2008 outburst while our high resolution spectroscopy was 
gathered during a quiescent phase. Our photometry conclusively determine that the outburst was due 
solely to the embedded Herbig Ae/Be member, supporting results from earlier works, and that the 
optically visible FU Ori component decreased slightly (~30%) in luminosity during the same period, 
consistent with previous works on the variability of FU Ori type systems. Further, our high-resolution 
ii'-band spectra definitively demonstrate that the 2.294 fim CO absorption feature seen in composite 
spectra of the system is due solely to the FU Ori component, while a prominent CO emission feature 
at the same wavelength, long suspected to be associated with the innermost regions of a circumstellar 
accretion disk, can b e assigned to the Herbig Ae/Be memb er. These findings are in contrast to 
previous analyses (e.g. lMalbet et aLll201Clt[Benistv et al.ll2010l ) of this complex system which assigned 
the CO emission to the FU Ori component. 

Subject headings: instrumentation: adaptive optics — stars: individual (Z CMa) — stars: pre-main 
sequence — binaries: close 



1. INTRODUCTION 

Z Canis Majoris (hereafter "Z CMa") was one of the 
original members of the set of Ae and B e stars with neb- 
ulosity first reported in iHerbid ()1960t ) , now called the 
Herbig Ae/Be type stars. This classification alludes not 
only to its early spectral type and emission lines, but also 
to brightness fiuc tuations and heavy nebulosity. At a dis- 
tance of 1150 pc (|Clarialll974D . Z CMa is pa rt of the CMa 
Tl a ssociation with a quoted age <1 Myr (jHerbst et al.l 
119781) . Consistent with its classification and presumed 
age, the Z CMa system shows strong P Cyg profiles 
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in fines of H^, H^ and Fe II (jCovino et al.lll98"^ . ra- 
dio continuum e mission, bipo l ar jet s as well as strong 
infrared excess ("P oetzel et al.l |1989() . Initial efforts to 
image the emitting so urce at high resolutio n revealed a 
~100 mas elongation (jLeinert fc Haaslll987[ ). suggesting 
a nearly edge-on dis k-like structure. Shortly thereafter, 
iKoresko et al.l (|1989D suggested that the extended emis- 
sion was probably not due to a disk, but rather due to 
a cooler, close companion retaining luminosity from its 
pre-main sequence contract ion. 

At about the same time, IHartmann et al.l (|1989D sug- 
gested the Z CMa system had many characteristics o f FU 
Ori-like systems (see e.g. IHartmann &: Kenvonlll996f) . in- 
cluding a bright accretion disk with a high rate of accre- 
tion, and signatures of powerful winds. Progress on rec- 
onciling the Herbig Ae/Be nature of the star with its now 
apparent FU Ori-like characteristics was achieved when 
evidence was uncovered through speckle imaging fo r 
two objects comprising the system fKores ko et al.|[l991l) . 
These m easurements were late r verified by visible direct 
imaging ()Thiebaut et al.iri995D revealing two clearly de- 
fined point sources. Such results provided support to the 
hypothesis that a significant fr action of FU Ori type sys - 
tems may be binary systems ([Bonnell fc Bastienlll992fl . 
which is indeed the case for the FU Ori system itself 
(|Wang et al.ll2004HReipurth fc Aspinll2004 Pueyo et al. 
2012). 

Synthesizing several decades of observations, a consis- 
tent picture of the Z CMa system has now emerged of a 
binary system comprised of a FU Ori-like object which 
dominates the optical emission, and a Herbig Ae/Be ob- 
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Fig. 1. — Three panels showing the Z CMa system shortly after its 2008 outbur st at 1.11, 1.44, and 1.71 /^m, obtained using the "Project 
1640" Integral Field Spectrograph at Palomar Observatory IHinklev et al.ll2011d) . The two members of this system are separated by ~1.5 
Palomar near-infrared diffraction widths. At longer wavelengths, the Northwest member, often referred to as a "Herbig Ae/Be" object, 
dominates the emission over the Southeast member, the "FU Orionis" object. 



ject, which dominates the in frared (jWhitnev et aL|[l993l : 
Ivan den Ancker et al.|[200^ . The terminology used for 
each component in the system in past hterature has been 
somewhat ambiguous given the constantly evolving phys- 
ical understanding of this remarkable system. As such, 
in this work we have chosen to abandon the "primary" 
and "secondary" nomenclature that pervades the litera- 
ture, and instead use the terms "Z CMa Southeast (SE)" 
to refer to the FU Ori component and "Z CMa North- 
west (NW)", referrin g the the Herbig Ae/Be component. 
iKoresko et alJ (|1991| ) show that the SE component dom- 
inates the emission of the NW component blueward of 
2/im, and vice versa for wavelengths longer than 2 /im. 

Support for the basic physical picture of the system 
has been verified repeatedly over the past ten years: i.e. 
the ide ntification that the SE component is a ~few Mq 
object (van den Ancker et al." '2004) responsible for the 
jet /outflow phenomena (Millan-Gabct & Monnicr 2002), 
while the NW component has considerable evidence 
supporting its Herbig Ae/Be nature incl uding emission 
lines and a compact, massive envelope (iWhitnev et al.l 
Il993| ). This syste m experienced a sigiiificant out- 
burst in 2008 (e.g. iGrankin fc Artemenkol l2009t ) last- 
ing '^l.S years, and was photometrically m onitored ex- 
tensi vely at X-ray to visible wavelengths (|Malbet et al.l 
20101). Indeed , usin g spectropolarimetry of the system. 



Szeifert et al.l (|2010i) claimed the outburst was associ- 



ated with the Herbig Ae/Be component of the system. 
Further, while dual-imaging polarimetry has typically 
been reserved for young circum stcllar disk sys tems (e.g. 
Qppenheimer et all 120081: iHinklev et al.ll2009f ). recently 



Canovas et al.l ( 20121) have used this technique to iden- 



tify the jets associated with each member of this system. 

Near-infrared spectroscopy of the Z CMa system 
has been scarce over th e past 30 years. In- 
deed, as lAnt oniucci et al. (2009) point out, only 
iLiliestroem fc O lofsson (1997) and v an den Ancker et al.l 
( 20041) have provided previous near-infrared spectra. In 
this paper, we present near-IR spectra and/or photom- 
etry in the JHKL bands for both components during 
the 2008 outburst, resolved using high angular reso- 



lution imaging (jOppenheimer fc HinklevI I2OO90 . as well 
as archival photome try covering the same bandpasses 
(|Koresko et al.lll991l) for comparison. Our adaptive op- 
tics ( "AO" ) measurements directly demons trate and pro- 
vide confirmation of previous work (e.g. iSzeifert et al.l 
EOIO; Malbct et al. 2010) that suggested the outburst is 
due to the Herbig Ac/Be component, although the exact 
mechanism of the brightening of this member remains 
unknown. Moreover, we present high-resolution K-hand 
spectra of each member of the binary obtained in 2006 
during the quiescent phase. Table [T] provides a list of our 
observations. 

2. OBSERVATIONS AND DATA PROCESSING 

Our high angular resolution observations of Z CMa 
were taken using Palomar/P1640 and Keck/NIRC2 AO 
during outburst. Our AO-assisted X-band spectroscopy 
was taken with Keck/NIRSPEC in quiescence. Table [T] 
provides a summary. 

2.1. JHKL-band Spectra & Photometry During 
Outburst 

We imaged the Z CMa system in its high state on 
UT 20 09 Marc h 17 using "Proj ect 1640" (jHinkle v et all 
l2011d :fOppcnh eimer et al]l2012h on the 200-in Hale Tele- 
scope at Palomar Observatory. Project 1640 is a coro- 
nagraph integrated with an integral field spectrograph 
(IFS). The IFS-I-Coro nagraph package is mounted on the 
Palomar AO system (|Dekanv et al.lll998l ). which in turn 
is mounted at the Cassegrain focus of the Hale Telescope. 
The cor onagraph is an Apodized-Pupil Lyot coron agraph 
(APLC; iSivaramakrishnan et al.ll25oiL ISoummed 12005) . 
having a 370 mas diameter {5.37X/D at i?-band) fo- 
cal plane mask. The IFS, is a microlens-based imaging 
spectrograph which can simultaneously obtain ^40,000 
spectra across our 3.8" x 3.8" field of view. Each mi- 
crolens subtends 19.2 mas on the sky and a dispers- 
ing prism provides a spectral resolution (A/AA)~45. 
The observing wavelengths spanned the J and H- 
bands (1.06/im - 1.76/im at the time of these observa- 
tions). Early examples of spectrophotometry a nd as- 
trometry from this project can be found in .Hinklev et al.l 
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(l2010l).lZimmerman et a ll (l20Tol) . iHinkley et all (l2011bi). 
[Roberts et all (|2012[ ). and lPuevo et alj (|2012bl ). While 



Fig. 2. — Left: Spectral Energy Distributions of each component of the Z CMa system in the quiescent state taken from lKoresko et al.l 
11199 II ). The NW (Herbig Ae/Be) component is shown with the red points while tlie photometry for the SE member (FU Ori) is shown 
in blue. Right: In the same color scheme, 1.06 — 1.76/.tm outburst spectra obtained from the Palomar Project 1640 spectrograph on 2009 
March 17 (open circles), as well as HKL Keck/NIRC2 points obtained on 2008 December 21 indicated by star-shaped points. Each solid 
curve represents two-component models comprised of the sum of two reddened blackbodies shown by dotted and dashed lines. The 2008 
outburst of this system is conclusively due to the order-of-magnitude brightness increase by the Herbig Ae/Be component. 

of the uncertainties in the J and 7?-band spectra. 

To avoid saturation on such a bright target, the out- 
burst data obtained with NIRC2 at Keck Observatory on 
UT 2008 December 21 utihzed the Jcont, ^^cont, i^cont, 
and PAH narrow filters with central wavelengths 1.213 
/im, 1.580 /im, 2.271 /xm, and 3.290 ^m, respectively. 
The flux values in these narrower filters were then scaled 
to corresponding broadband JHKL wavelengths using 
the ratio of bandpasses between the two sets. The Z 
CMa system was observed using three-point dithering, 
and airmass corrected using the calibration stars listed 
in Table □ 

The left panel of Figure [2] shows the photometry of 
each component in the quiescent state collected from 
'Kores ko et al.l (|199H) using speckle interferometric obser- 
vations. The right panel of Figure [5] shows our J and H- 
band i?^45 spectra from Palomar, as well as our HKL- 
band photometry obtained at Keck for both members of 
the system taken from the 2008-2009 outburst observa- 
tions. 



an IFS clearly aids in spectral characterizat ion, it can 
also i mprove sensitivity (Crepp et al. 2011; P uevo et al.l 
l2012ail through the suppression of quasi-static speckle 
noise in the direct image which w ill limit high contrast 
observations (|Hinklev et al.|[2"007[ ). This system is also 
equipped with aperture maski ng interferometry capa bil- 
ities (e.g. lHinklev et al.ll2011al : lKraus fc Irelandi[2012l ) for 
the characterization of systems with very small angular 
separations, although this technique was not employed 
in this study. 

Three individual wavelength channels are shown in 
Figure [T] To fully resolve both members of the system, 
the target was observed ~1 arcsec away from the coron- 
agraphic mask. 

Point Spread Function ("PSF") fitting photometry was 
performed in each of the 23 wavelength channels using 
a calibration star, HD 112196 (F8V, V=7.01), obtained 
on the same night with similar signal-to-noise and ob- 
served at similar airmass. The 100 mas angular sepa- 
ration of the two Z CMa components (Figure [IJ is only 
slightly larger than the 70 mas _ff-band diffraction limit 
of the Palomar Hale Telescope. Given this fact, fitting 
the empirical calibration PSF star to each member was 
performed simultaneously. Each channel of the Z CMa 
images were oversampled by a factor of five, and the 
calibration PSF was fit to each Z CMa component to 
determine peak brightness and centroid positions. Once 
these two best-fit PSFs were determined, each was subse- 
quently subtracted from their corresponding Z CMa com- 
ponent to determine the post-subtraction residual fiux. 
This residual fiux is the dominant term in the calculation 



2.2. High Resolution K-band spectroscopy During 
Quiescence 

Spatially resolved if-band spectra of the two Z CMa 
components were obtained with the NIRSPEC spectro- 
graph mounted behind the ad aptive optics system on 
the W. M. Keck II telescope (jWizinowich et all l2000( ) 
on UT 2006 December 17. The if-band setup provided 
wavelength coverage across the CO bandhead, the Br 7 
line, a He I line, and a Mg I and Al I metallic line re- 
gion, but the spectral coverage is not continuous. The 
i?=A/AA'-^25,000 observations were obtained at two po- 
sition angles: 60 degrees (along the jet axis) and 120 
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TABLE 1 

Table of Observations 



Object Star Type 



Date 



Z CMa State Observatory & Instrument Wavelengths (^tm) 



Mode 



Z CMa 
S Mon 
Z CMa 
HD 53455 
HD 75898 
Z CMa 



Target 

Calibrator 

Target 

Calibrator 

Calibrator 

Target 



2006 Dec 17 
2006 Dec 17 
2008 Dec 21 
2008 Dec 21 

2008 Dec 21 

2009 Mar 17 



Quiescence 
Outburst 



Keck 
Keck 
Keck 
Keck 
Keck 



NIRSPEC 

NIRSPEC 

NIRC2 

NIRC2 

NIRC2 



HD 112196 Calibrator 2009 Mar 17 



Outburst Palomar: Project 1640 
Palomar: Project 1640 



2.10-2.13, 2.29-2.32 
2.10-2.13, 2.29-2.32 

1.57-3.32 

1.57-3.32 

1.57-3.32 

1.06-1.76 

1.06-1.76 



Spectroscopy 
Spectroscopy 
Imaging 
Imaging 
Imaging 
Spectro-photometry 
Spectro-photometry 



degrees (along the projected semi-major axis of the bi- 
nary). The slit was 2.26" long and data were taken using 
two different widths: 0.068" and 0.027". MCDS readout 
mode was used with 30s and 60s integration times with 
five co-adds. Because of the complexity of the field, in 
addition to dithered pairs of spectra taken in an ABBA 
nod pattern, off-target sky images were also obtained. 
Telluric correction was achieved using the rapidly rotat- 
ing early type star S Mon and further calibration was 
derived from internal arc line lamps and fiat field lamps. 
Figure [3] shows the flux from two echelle orders from 
NIRSPEC, spanning 2.10 - 2.13^j.m and 2.29 - 2.32/m. 

3. RESULTS 

3.1. Confirmation of the Outbursting Component 
The quiesc ent flux measurements taken from 



iKoresko et alJ (|1991l ) and shown in Figure [2] reveal 
that shortward of ~2 /im, the overall system flux is 
dominated by the SE component (blue curve in Fig- 
ure [5]) . Longward of ^2 /im, the NW component (the 
Herbig Ae/Be star shown by the red points) dominates. 
Further, the NW component has increased in brightness 
during outburst by an order-of-magnitude over the 
quiescent values. At the same time, we find a '~30% 
flux decline between the quiescent and outburst states 
for the SE component (the FU Ori star), shown in blue. 
The high angular resolution afforded by the Keck and 
Palomar AO systems allows us to confirm the claim 
(iMalbet et all 120101 : iBenistv et all 120101 : iSzeifert et all 
120101 ) that the NW component is responsible for the 
ensemble system's increase in brightness. 

3.2. Two-component Model Fits to Data 

Using the photometry obtained at Keck and Palomar, 
as well as the work of Koresko ct al. (1991), we are able 
to satisfactorily fit two-component, reddened, blackbody 
models to the quiescent and outburst Spectral Energy 
Distributions for each member of the Z CMa system 
shown in Figure [2l For the NW component, our model 
that is physically plausible for such a luminous star is 
made up of two blackbodies of temperatures 8500 and 
HOOK behind Av=10. For the SE member, we find 
a physical match using components of 5500 and 900K 
behind Av=5.2. The reddening uses a o ptical total-to - 
selective extinction ratio, i?^=3.1, and the lMathii ()1990f ) 
extinction curve. Although we have chosen to refrain 
from assigning a great deal of physical meaning to this 
two-component model, we refer to these components as 
"photosphere" and "disk" . In each of the quiescent and 
outburst phases, the disk blackbody temperatures, nor- 
malizations, and reddening values are identical. At the 
same time, the normalization of the fit to the photo- 



spheric component of the NW member has increased by 
an order of magnitude from quiescence to outburst, while 
for the SE component, it has been reduced by 30% be- 
tween the two stages. The inherent variability of the Z 
CMa system is reflected in these changing photosphere 
normalizations. 

Figure [2] shows that in addition to the order-of- 
magnitude brightness increase of the NW component, 
the peak luminosity of this member has shifted blue- 
wards, indicating either a greater luminosity of the pho- 
tosphcric component, or a significant decrease in the Ay 
value. However, decreasing the Ay value immediately 
requires unphysical values for the photospheric temper- 
ature. Specifically, Ay— 2 for the NW component im- 
plies an unphysical best flt photosphere temperature of 
'^1300K. The 5500K photospheric temperature for the 
SE member agrees well with typical FU Ori objects that 
show CO absorption fe atures (See Section 13. 3|) as de- 
scribed, for example, in lCalvet et al.l ()1991[ ). Meanwhile 
the disk temperature for the NW membe r agrees well 
with values quoted in ISzeifert et al.l (|2010( ). Lastly, this 
~30% decrease in flux for the SE component over '^20 
ye ars is remarkably similar to the optical decay reported 
in lClarke et al.l ()2005i) for the archetypal system FU Ori. 

3.3. CO Features Observed in Quiescence 

The NIRSPEC spectra shown in Figure [3| capture the 
prominent He I emission line at 2.113 /im as well as the 
CO bandhead at 2.294 /im. CO emission in this wave- 
length range occurs in ^- ^20% of luminous young stellar 
systems (e.g. lCarij|1989[ ). and probes the i nnermost re- 
gions of the circumstellar accretion disks ()Naiita et alJ 
l2008h . 

Our Keck high angular resolution observations al- 
low us to determine that the SE component (the FU 
Ori star) possesses a significant CO bandhead absorp- 
tion at 2.294 /Ltm. Several FU Ori objects such as FU 
Ori itself, V1057 Cyg, and V1515 Cyg, show signifi- 
cant absorption features due to this CO band head (e.g. 
[Reipurth & Aspin 1997). It has been noted that Z CMa 
possesses rather weak CO absorption at 2.294 /im rela- 
tive to other FU Ori stars, attributed to mixing of the 
signals between th e two components due to in adequate 
angular resolution (jHartmann fc Kenvo nl 119961 ) . We can 
check this hypothesis with our high angular resolution 
observations. Our clear detection of the absorption band 
head at 2.294 /im in the SE component confirms the 
claim of Ivan den Ancker et al.l (j2004l) that the absorp- 
tion is most likely arising fro m this component. While 
iLiliestroem fc OlofssonI (|1997[ ) present evidence for the 
CO band head using medium resolution K-hand spec- 
tra (see especially their Figure 1), the spectral resolution 
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Fig. 3. — Keck K-hand spectra from the NIRSPEC spectrograph 
for both the SE and NW components of Z CMa obtained while the 
system was in quiescence in 2006 December 17. The Herbig Ae/Be 
component shows a prominent CO emission feature located at 2.293 
/jm, while the FU Ori component shows a significant absorption at 
this wavelength. 

of their observations was not sufficient to ascribe these 
CO bandhead features to one component or the other. 
Indeed, these authors attempt a decomposition of the 
spectrum by choosing a hnear analytic function for the 
NW component, which forces deviations from this linear 
form to be due to the SE component. 

At the same time, Figure [3] shows that the NW com- 
panion (the Herbig Ae/Be) s hows an emis s ion fe ature 
at the band head location. iMalbet et al.l (|201C1( ) and 
iBenistv et all ()2010D note the presence of this CO over- 
tone line in emission during outburst in their spatially 
unresolved observations of the Z CMa system, but as- 
sume that this arises from the SE component (the FU 
Ori member). Our high spectral resolution observations 
clearly demo nstrate the CO emiss ion arises from the NW 
component. iCalvet et al.l ()199lD discuss at length the 
conditions under which an object may exhibit CO emis- 
sion during an outburst phase through the increased ir- 
radiation of the atmosphere of the circumstellar disk. 

4. CONCLUSIONS 

We present high angular resolution near-infrared ob- 
servations of the Z Canis Majoris system during its 2008 



outburst as well as prior to this during a quiescent phase. 
The JHKL-hand outburst photometry conclusively de- 
termines 1) that the brightening is due solely to the em- 
bedded Herbig Ae/Be member, confirming results from 
earlier works, and 2) that the optically visible FU Ori 
component has actually experienced slightly declining 
brightness between the quiescent and outburst stages. 
Thus, this substantial brightness increase of the Herbig 
Ae/Be component relative to its quiescent state is re- 
sponsible for the continuum brightness doubling of the 
ensemble system. Further, the Keck high-resolution K- 
band spectra obtained during a quiescent phase allow 
us to conclusively determine that the 2.294 iim CO ab- 
sorption feature seen in composite spectra of the sys- 
tem is due solely to the FU Ori component. In addition 
however, these X-band observations show a prominent 
CO emission feature in the Herbig Ae/Be member which 
likely dilutes the strength of the corresponding absorp- 
tion feature in the SE component when the images of 
each component are blended with lower angular resolu- 
tion observations. This result stands in contrast to as- 
sumptions in oth er works a bout the source of the CO 
emission (e.g. IBenistv et al.lfeoiO.: .Malbet et al..,201Q) . 
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